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Abstract 


Fuzzy Logic Controllers (FLCs) using Genetic Algorithm have been designed for four 
different dynamical systems. The parameters of fuzzy logic controller have been deter- 
mined optimally by using genetic algorithm with multi objective cost function. 

1. Poor water level control of U tube steam generator (UTSG) in nuclear power 
plant may lead to frequent nuclear shutdown. UTSG of pressurized water reactor 
(PWR) has a time varying dynamics. FLC has been designed for this system and 
its performance is compared with LQR controller. It is found that its performance 
is as comparable with LQR controller. 

2. Control of a U tube steam generator is possible if a proper mathematical model 
can be derived. A U tube steam generator has been modeled for analyzing its 
behavior. The PD, PI and fuzzy logic controllers have been implemented for 
controlling the water level. 

3. FLC controller has been implemented on a single link manipulator and its perfor- 
mance has been tested with linearized and PD controllers. It is found that FLC 
performs better than PD Controller and is well compared with linearized control 
schemes. 

4. FLC has been used to balance a reaction wheel pendulum, which is an under 
actuated system. 
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Chapter 1 


Introduction 

1.1 Introduction 

All natural systems are inherently nonlinear. The properties of nonlinear systems are 
quite different from those of linear systems. For example, nonlinear systems exhibit 
properties like limit cycles, chaos, bifurcation etc. Moreover its stability depends not 
only on the initial condition but also on the control input. Linear system theory is very 
well developed and numerous tools are available for its analysis. In nonlinear systems, 
the tools available are limited and they are problem specific. Various linear system 
techniques have been used to analyze and control the linearized model of nonlinear 
plants. Recently, a number of other tools belonging to class of “learning systems” have 
been developed for designing high performance controllers. This later class of tools 
include Fuzzy logic, neural networks, evolutionary algorithms etc. 

This thesis consists of two parts. In first part, we discuss about the modeling and 
control of a nuclear boiler and in the other, we discuss about various controllers for 
nonlinear mechanical systems. 
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1.2 Level Control of Steam Generators 

Nuclear power is identified as one of the major electric power source in future. These 
power plants are very sensitive to any disturbance (mechanical and electrical) due to 
various technical and safety reasons. U tube steam generator is a major part of nuclear 
power plant. It forms a vital link between the heat source (the reactor) and the heat sink 
(the turbo generator system). The prediction of this behavior requires the knowledge 
of non-linear dynamical models. 

Controlling the water level in a large PWR plant is of major importance to ensure 
sufficient cooling of reactor, good performance of steam separators and dryers. Simulta- 
neously it is required to eliminate all risks of hydrodynamic instability [1]. As a result, 
a good control system is of paramount importance for the proper functioning of the 
nuclear plant. 

Since the actual plant is nonlinear and complex, most of the time, linearized model 
is used for designing controllers. For linearized plants, controllers like LQR and state- 
feedback control have been used extensively in the literature [2] [3]. Keeping the com- 
plex nature of boiler dynamics, a number of techniques based on soft computing ap- 
proaches like “Fuzzy Logic” have been developed for designing controllers. [4] [5] . Fuzzy 
logic controllers have been shown to improve the performance of existing conventional 
controllers. 


1.3 Control of nonlinear mechanical systems 

Nonlinear dynamical system like single link manipulator, reaction wheel pendulum are 
used as a bench mark problem for control system. Some systems are under actuated 
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which are more difficult to control. For nonlinear plant, a number of different controllers 
have been used in the literature [6] [7] [8]. Different non linear control schemes like 
adaptive control, fuzzy logic control, and sliding mode control have been described in 
the literature. In the present work fuzzy logic control is used to control three non linear 
control systems. 

1.4 Contribution 

The performance of FLC has been compared with those of conventional controllers for 
different class of nonlinear systems. The nonlinear systems studied here, include single 
link manipulator, reaction wheel pendulum and steam generator. The fuzzy parameters 
have been optimized using genetic algorithm. 

Detailed modeling of a U tube steam generator has been carried out based on the 
technical data provided by BARC. 

A fuzzy logic controller has been designed for water level control of U tube steam 
generator in a nuclear power plant and its performance is compared with those of 
conventional controllers. 

All the simulation works have been carried out in C++ language. 

1.5 Thesis Organization 

The thesis is organized as follows. 

Chapter 2 describes briefly about Fuzzy Logic and Genetic Algorithm. 

In Chapter 3, LQR and Fuzzy Logic controllers are used for level control of U tube 
steam generator. Here we use a linearized model of the actual plant. 
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The modeling of Nuclear steam generator has been carried out in Chapter-4. Three 
different controllers have been implemented on this model. 

In Chapter 5, performance of different controllers are compared for a single link 
manipulator problem. 

The balancing problem of a reaction wheel pendulum is analyzed in Chapter 6. 
Performance of different controllers are compared for balancing task. 

The thesis is summarized and its future scope is discussed in Chapter 7. 

The overview of different controllers and formulae related to Boiler dynamics have 
been provided in the Appendix. 
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Chapter 2 


Introduction to Fuzzy Logic and 
Genetic Algorithm 

2.1 Introduction to Fuzzy Logic 

In the real world, information is often ambiguous or imprecise. When we state that it is 
warm today, the context is necessary to approximate the temperature. A warm day in 
January may be -5 degrees Celsius, but a warm day in August my be 35 degrees. After 
a long spell of frigid days, we may call a milder but still chilly day relatively warm. 
Human reasoning filters and interprets information in order to arrive at conclusions or to 
dismiss it as inconclusive. Although machines cannot yet handle imprecise information 
in the same -ways that humans do, computer programs with fuzzy logic are becoming 
quite useful when the sheer volume of tasks defies human analysis and action. 

An organized method for dealing with imprecise data is called fuzzy logic. The data 
are considered as fuzzy sets. For traditional sets, there is no other case than true or 
false. Fuzzy sets allow partial membership. Fuzzy Logic is basically a multi valued 
logic that allows intermediate values to be defined between conventional evaluations 
like yes or no, true or false, black or white etc. Notions like rather warm or pretty cold 
can be formulated mathematically and processed with the computer. In this way an 
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attempt is made to apply a more human-like way of thinking in the programming in 
the computers. 

2.2 Fuzzy Logic - Its Origin 

The concept of fuzzy logic (FLC) was conceived by Lotfi Zadeh, a professor at the 
University of California at Berkley, and presented not as a control methodology, but 
as a way of processing data by allowing partial set membership rather than crisp set 
membership or non-membership [4]. Fuzzy logic was not applied to control systems 
before 70’s due to insufficient micro computer capability prior to that time. Professor 
Zadeh reasoned that people do not require precise, numerical information input, and 
yet they are capable of highly adaptive control. 

2.3 How Fuzzy Logic is Different from Conventional 
Control Methods? 

FLC incorporates a simple, rule-based IF X AND Y THEN Z approach for solving 
control problem rather than attempting to model a system mathematically. The FLC 
model is empirically based, relying on an operator’s experience rather than their tech- 
nical understanding of the system. For example, rather than dealing with tempera- 
ture control in terms such as “TEMP =500 °F”, ’’TEMP < 1000 °F ”, or ”210 °F 
<TEMP< 220 °F ” , terms like “IF (process is too cool) AND (process is getting colder) 
THEN (add heat to the process)”. FLC is capable of mimicking this type of behavior 
but at very high rate. 
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2.4 How Does Fuzzy Logic controller Work? 


Fuzzy Logic controls requires some numerical parameters in order to operate such as 
what is considered significant error and significant rate-of-change-of-error, but exact 
values of these numbers are usually not critical unless very responsive performance is 
required. 

2.5 Here are some guide lines by which a Fuzzy logic 
controller is designed 

1. Define the control objectives and criteria: What is the control objectives? What 
are the parameters of the system to be controlled? What kind of response is 
needed? What are the possible (probable) system failure modes? 

2. Determine the input and output relationship and choose a minimum number of 
variables for input to the FLC engine (typically error and rate of change of error 
for fuzzy PD Controller). 

3. Using the rule-based structure of FLC, break the control problem down into a 
series of IF X AND Y THEN Z rules that define the desired system output 
response for a given system input conditions. Although it is possible to use a 
single, instantaneous error parameter without knowing its rate of change, this 
cripples the system’s ability to minimize overshoot for a step inputs. 

4. Create FLC membership functions that define the meaning (values) of Input/Output 
terms used in the rules. 
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5. Test the system, evaluate the results, tune the rules and membership functions , 
and continuously simulated until satisfactory results are obtained. 

2.6 Architecture of Fuzzy Logic Controller (FLC) 

A general architecture of fuzzy logic controller is shown in figure 2.1. Various subsystems 
of FLC are described in the following sections. 



(Fuzzy Controller Architecture) 

Figure 2.1: Fuzzy Architecture 


2.6.1 Fuzzification 

As shown in the figure 2.1, fuzzification is a process in which by the crisp inputs 
are mapped into corresponding linguistic (fuzzy) variables. These linguistic variables 
are used for making decisions. Each linguistic variable or Fuzzy set has a certain 
membership function. The membership function assigns a certain degree of membership 
or weight for each of the linguistic variables. The membership value is nothing but the 
truth of belongingness of a certain input to that fuzzy set or linguistic variable. 


8 






2.6.2 Designing the rule base 


IF-THEN Rules 

Fuzzy sets and fuzzy operators are the subjects and verbs of fuzzy logic. These IF- 
THEN rule statements are used to formulate the conditional statements that comprise 
fuzzy logic. 

A single fuzzy IF-THEN rule assumes the form if x is “A” then y is “B” . Where “A” 
and “B” are linguistic values defined by fuzzy sets on the range (universes of discourse) 
x and y respectively. The if-part of the rule “x is A” is called the antecedent or premise, 
while the then part of the rule y is “B” is called the consequent or conclusion. 

A double fuzzy IF-THEN rule base is if x is “A” , y is “B” than z is “C” . A sample 
81 fuzzy rule base is shown in the Table 2.1. First row of table indicate fuzzy sets for 
first variable (A). First column of the table indicate fuzzy sets for second variable ( B ). 
Other elements of the table indicate rule for the generating control variable. 

Where nvl, nl, nm, ns ze, ps, pm, pi, pvl are the fuzzy sets for A named as negative 
very large, negative large, negative medium, negative small, zero, positive small, positive 
medium, positive large, positive very large respectively, dnvl, dnl, dnm, dns dze, dps, 
dpm, dpi, dpvl are the fuzzy sets for B named as negative very large, negative large, 
negative medium, negative small, zero error, positive small, positive medium, positive 
large, positive very large respectively, unvl, uni etc are the rules for generating control 
input. Each of these variable are assigned with a specified value. All these variables 
are assigned by a specific value. 
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rules 

nvl 

nl 

nm 

ns 

ze 

ps 

pm 

pi 

pvl 

dnvl 

unvl 

unvl 

unvl 

unl 

unl 

unm 

unm 

uns 

uze 

dnl 

unvl 

unvl 

unl 

unl 

unm 

unm 

uns 

uze 

ups 

dnm 

unvl 

uni 

unl 

unm 

unm 


uze 

ups 

upm 


uni 

unl 

unm 

unm 

uns 


ups 

upm 

upm 

EH 

uni 

unm 

unm 

uns 

uze 

ups 

upm 

upm 

upl 

dps 

unm 

unm 

uns 

uze 

ups 

upm 

upm 

upl 

upl 

dpm 

unm 

uns 

uze 

ups 

upm 

upm 

upl 

upl 

upvl 

dpi 

uns 

uze 

ups 

upm 

upm 

upl 

upl 

upvl 

upvl 

dpvl 

uze 

ups 

upm 

upm 

upl 

upl 

upvl 

upvl 

upvl 


Table 2.1: FAM Table for Rule Base 


2.6.3 Fuzzy Inference Mechanism 

Fuzzy inference is the process of formulating the mapping from a given input to an 
output using fuzzy logic. The mapping then provides basis from which decisions can be 
made, or patterns discerned. Now control decisions can be made based on the linguistic 
variables. It performs two functions. Firstly assigns degree of membership or weight to 
the variables and secondly determines the output decisions or output linguistic variables 
based on input linguistic variables. 


2.6.4 Defuzzification Model 

T his converts the fuzzy set of output variables to crisp numerical outputs. These crisp 
values are the final control inputs to the system to be controlled. 


2.7 FLC Design 


The design process of an FLC may be divided into the following steps described below. 
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2.7.1 Selection of Control Variables 


The selection of control variables (controller inputs and outputs) depends on the nature 
of the controlled system and the desired outputs. It is more common to use the output 
error (e) and the rate or derivative of the output (e) as controller inputs. 

2.7.2 Membership Function 

Each of the FLC input signals, output signals and fuzzy variables, have the real line 
“R” as the universe of discourse. In practice, the universe of discourse is restricted to a 
comparatively small interval [Xminj, Xmaxj]. The number of fuzzy sets for each fuzzy 
variable varies according to the application. A common and reasonable number of vari- 
ables are odd numbers. Increasing the number of fuzzy sets results in a corresponding 
increase in the number of rules. A membership function is assigned to each fuzzy set 
[9]. The membership function maps the crisp set into fuzzy set. A set of membership 
functions defined for nine linguistic variables NVL, Nl, NM, NS, Z, PS, PM, PL, PVL 
which stand for negative very large, negative large, negative medium, negative small, 
zero, positive small, positive medium, positive large, positive very large respectively. 

Membership function can be of a variety of shapes, the most popular structures are 
defined as triangular, trapezoidal, or a bell shaped. For simplicity, it is assumed that 
the membership functions are symmetrical and each one overlaps with the adjacent 
functions by 50%. In practice, the membership functions are normalized in the interval 
[-L, L], which is symmetrical around zero. Thus control signal amplitudes (fuzzy vari- 
ables) are expressed in terms of controller parameters (gain). A triangular membership 
function is given in the figure 2.2 
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-20 -15 -10 -5 0 5 10 15 2 0 

Input space for x=[-20, 20] 


Figure 2.2: Triangular Membership Function for 9 Regions 
2.7.3 Defuzzification Strategy 

Defuzzification is a process of converting the FLC inferred control actions from fuzzy 
values to crisp values. This process depends on the output fuzzy set, which is gener- 
ated from the fired rules . The output fuzzy correlation is formed by either minimum 
encoding or the product encoding. 

Commonly defuzzification methods have been discussed below. 

Centroid Method 

This procedure also known as center of area or center of gravity method is the most 
prevalent and physically appealing of all the defuzzification methods. It is given by the 
algebraic expression 

EIU 
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Weighted average method 


This method is only valid for symmetrical output membership functions. It is formed by 
weighting each membership function by respective membership value. This is written 
in mathematical expression as. 

E n 

1 V— \7l 

2_/i=l fa 

where 

n=number of rules. 

/^membership for i t h rule. 
cu=value of that particular rule. 

T— Control variable. 

Max membership principle 

T his method also known as height method. This scheme is limited to peaked output 
function. This method is given by the algebraic expression 
fair) > m{u) 


2.8 Summary 

FLC offers several unique features that make it good choice for many control problems. 

• FLC was conceived as a better method for sorting and handling data. It has 
proved to be a excellent choice for many control system applications since it 
mimics human control logic. It can be built into anything from small, hand made 
products to large computerized process control systems. 
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• Fuzzy based control (FLC) has become highly competitive due to its better per- 
formance, high reliability, robustness,' low power consumption and cheapness. 

• The thinking process involved in fuzzy realm is not complex. It is simple, elegant 
and easy to apply. 

• In this context, FLC is one of the methodologies for solving control system prob- 
lem. It lends itself for implementation in systems ranging from simple, small, 
embedded micro controllers to large, networked, multi-channel PC or workstation- 
based data acquisition and control systems. It can be implemented in hardware, 
software, or a combination of both. 

• FLC provides a simple way to arrive at a definite conclusion based upon vague, 
ambiguous, imprecise, noisy, or missing input informations. FLC approach to 
control problems mimics how a person would make decisions at a much faster 
rate. 

• It is inherently robust since it does not require precise, noise-free inputs. The 
output control is a smooth control function despite a wide range of input varia- 
tions. 

• Since the FLC processes user defined rules governing the target control system, it 
can be modified and tweaked easily to improve or alter the system performance 
drastically. New sensors can easily be incorporated into the system simply by 
generating appropriate governing rules. 

• FLC is not limited to a few inputs or few outputs. This allows the sensors to 
be inexpensive and imprecise. It keeps the overall system cost effective and low 


14 



complexity. 


• Because of the rule-based operation, any reasonable number of inputs can be 
processed and numerous outputs can be generated. Although defining the rule 
becomes complex if too many inputs and outputs are chosen for a single imple- 
mentation. It would be better to break the control system into smaller chunks 
and use several smaller FLC that can be distribute on the system. 

• FLC can be used to control nonlinear systems that would be difficult or impossible 
to model mathematically. This opens doors for control systems designers that 
would normally be deemed infeasible for automation. 

• FLC provides a different approach to control complex systems. This method 
focuses on what the system should do rather than trying to model how it works. 
One can concentrate on solving the problem rather than trying to model the 
system mathematically, if that is even possible. 

• On the other hand the fuzzy approach requires a sufficient expert knowledge for 
the formulation of the rule base, the combination of the sets and the defuzzifica- 
tion. 


2.9 Genetic Algorithm (GA) 

It is adaptive heuristic search algorithm premised on the evolutionary ideas of natural 
selection and genetic. The basic concept of GA is to simulate processes in natural 
system necessary for evolution, specifically those that follow the principles first laid 
down by Charles Darwin of “survival of the fittest” . As such they represent an intelligent 
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exploitation of a random search within a defined search space to solve a problem. 

Genetic algorithm is a search and optimization method developed by mimicking 
the evolutionary principles and chromosomal processing in genetics. It is a part of 
evolutionary computing which is a rapidly growing area of artificial intelligence. GA 
uses a chromosomal representation which required the solution to be coded in a finite 
length string. 

It follows the following steps for minimizing the cost function by optimizing some 
parameters. 

1. Representation. 

2. Reproduction. 

3. Crossover. 

4. Mutation. 

The flow chart of genetic algorithm is shown in figure 2.3. 

2.10 Genetic Algorithm for Single Objective Opti- 
mization 

Genetic algorithm identifies the parameters and its limit which are to be optimized for a 
given cost function (performance index). An example for a single objective optimization 
is illustrated below. 

Within the range of selection around 200 children are generated randomly by taking 
help of binary conversion. Out of these around 20 better children are selected for 
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Figure 2.3: Basic Flow Chart of Ganetic Algorithm 

minimizing objective function. By permutation and combination again 200 different 
population are generated for next generation by the process of crossover. Mutation is 
performed with less density than crossover for finding out new search. The number 
of generation is found till all the children have same chromosome (same value). Then 
it is found that no further development of the performance index is possible. The 
parameters at this performance index are chosen to be optimized value. 

2.10.1 Advantages 

It uses the natural selection process. For every generation, it takes better values of the 
parameter for next generation. So it is completely focus based search. It uses crossover 
and mutation for searching best performance. 
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2.10.2 Limitation 


It can search up to variation of 3 parameters. Some times it shows local minima when 
limit of the parameters are large. This can not be used for more than one objective 
function where most of the problems are multi objective function. 

2.11 Nsga2.code: 

It is used for solving a multi objective optimization. It is a code available in the web 
site of “Dr. K. Deb, IIT Kanpur” . In my Problem it has been tried to set the parameter 
such that it minimize the error and the control input throughout the simulation. 


18 



Chapter 3 

Level Controller for a U Tube 
Steam Generator 


3.1 Introduction 

Poor control of steam generator water level of a nuclear power plant may lead to frequent 
nuclear reactor shut downs. This shutdown is more common at low power operation. 
There is a need to investigate this problem so that a suitable controller for water level 
regulation can be designed. This work is concerned to study and design of a suitable 
controller for a U Tube Steam Generator (XJTSG) of a Pressurized Water Reactor 
(PWR). It has a time varying dynamics. 

One of the objectives for a steam generator in a nuclear power plant is to control the 
water level at a desired value by regulating the feed water flow rate. The conventional 
feed water control scheme cannot provide satisfactory performance within the required 
wide operating range of 0 to 100% of the specified load. A large proportion of the reactor 
shutdown causes a severe economic loss. It has been mainly caused by ineffective feed 
water control. Therefore a better control scheme is required. Past experiences on PWR 
plant operation, had indicate that frequent reactor trips are due to poor control of water 
level in steam generator. During low power operation, the level control is complicated 
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by the thermal reverse effects known as ” shrink and swell”. Due to the destabilizing 
vapor content in the tube bundle region, the water level measured in the down comer 
temporarily reacts in a reverse manner to water inventory change. Increased feed water 
flow adds mass to the SG, which would be expected to increase the measured down 
comer water level. It does not increase at high power. But at low power, the cold feed 
water addition can cause a decrease in the vapor content of tube bundle. A shift in the 
liquid from the down comer to tube bundle causes a temporary decrease in water level 
(shrink). Similarly, a decrease in feed- water flow can cause a temporary increase in 
water level (swell). These reverse effects are confusing for either manual or automatic 
operation. The only true indication of water inventory change is the flow mismatch 
between steam and feed-water. 

The water level of the steam generator must not be allowed to rise too high, in order 
to prevent the excessive moisture carryover and the pressure buildup of the containment 
due to discontinuity of secondary side flow loop. Also, the low water level should be 
prevented in order to avoid the uncovery of the U-tubes in the secondary side. Therefore, 
the control of the steam generator water level is important to determine power plant 
dynamics due to change in operating load. 

For the given system FLC and LQR controller are designed and simulated. There 
are some reports on automatic operation of the steam generator water level from low 
power to high power range [2] . 
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3.2 Model Development 

A steam generator model is a complicated nonlinear dynamical system. The controller 
design and its performance in actual plant strongly depend on the accuracy of the 
mathematical model. A highly accurate model of a steam generator is complex and 
nonlinear. So it is desired that the model should be simple and at the same time 
relatively accurate for describing the process of dynamics of U Tube Steam Generator 
(UTSG). 

The difficulty in designing the effective water level control for the steam generator 
arises from a number of factors which are given below. 

(a) The inverse response of the plant at low operating power due to the “swell and 
shrink” . 

(b) Variation of the plant dynamics with the operating power. 

(c) Unreliable flow measurements at low power which preclude effective use of feed- 
forward control. 

We have taken a model which is widely used by many researcher for the control purpose. 
This is a fourth order model whose parameters are dependent on reactor power level. 
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3.3 Dynamic Model of PWR for Level Measure- 
ment 

Transfer function of the water level control system of a UTSG is shown in the figure 
3.1. It can be written in mathematical form as, 


,,, „ Gi(q w (s) - q v (s)) g 2 (g^ (g) ~ 9v( s )) , G 3 sq w (s) 

r (a) 1 1 _ „ ~ -2 


1 + r 2 S 


t ,~ 2 + 4tt 2 T- 2 + 2r{s + s 2 


( 3 . 1 ) 


where, 

Y(s) : Water level in meter. 

q w : Feed water flow rate in kg/sec. 

q v : Steam flow rate in kg/sec. 

Ti,T 2 : Damping time constants in second. 
T : Oscillating period in second. 
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G\/ s is the mass capacity effect of the UTSG. It integrates the flow error Q w (s)—q v (s) 
to calculate the change in water level. This term accounts for the level change due to 
feed water inlet to steam generator and the steam outlet from it. This quantity Gi/s is 
the actual water capacity which critically affects the removal capability of the primary 
heat. Gi is a positive constant and does not depend on load. G 2 /( 1 + r 2 s) is the 
thermal negative effect caused by ” swell and shrink”. Since these phenomena exhibit 
exponential responses for step changes of the feed water flow rate and the steam flow 
rate corresponding dynamics are described by a first order equation. (? 2 is a positive 
constant and is dependent on load. As load increases (? 2 decreases. The third term in 
equation (3.1) is the mechanical oscillation effect caused by the inflow of the feed- water 
to the U tube steam generator. This is a mechanical oscillation term due to momentum 
of the water in the down comer. All the water removed from the steam is returned 
to the down comer and is recirculated. The recirculating water has large momentum 
acting against relatively small change in flow rate. When the feed water flow rate is 
suddenly decreased, the water level in the down comer falls initially and then begins 
to oscillate. This is due to the momentum of the water in the down comer keeping 
the recirculating flow going down initially and then slowing down. The mechanical 
oscillation disappears completely after a small multiple of the damping time constant. 
The variable Gz is a positive constant and depends upon the load. 

The steam generator dynamics has been divided into four linearized regions with 
respect to operating power level. It is assumed that the dynamics vary linearly over 
these regions. These variations of the plant parameters with respect to power level are 
presented in the table (3.1). The actual plant parameters may vary differently, so -we 
study the performance of the controller under the situations when the parameter drifts 
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q v (kg/second) 

57.4 

180.8 

381.7 

660 

1435 

P(%power in Mw) 

5 

15 

30 

50 

100 

G x 

0.058 

0.058 

0.058 

0.058 

0.058 

G 2 

9.63 

4.46 

1.83 

1.05 

0.47 

Gz 

0.181 

0.226 

0.31 

0.215 

0.105 

n 

41.9 

26.3 

43.4 

34.8 

28.6 

T 2 

48.4 

21.5 

4.5 

3.6 

3.4 

T 

119.6 

60.5 

17.1 

14.2 

11.7 


Table 3.1: Dynamic parameters with respect to operating power 


from what is projected by linear interpolation. 


3.4 Controller Design for Water Level Control 

Systematic approach for linear quadratic regulator (LQR) and fuzzy logic methods have 
been used to derive the control law where some objectives functions are minimized to 
derive an optimal controller. The main objective of the controller is to maintain the 
water level in the steam generator under various operating levels. We also show the 
effect of parameter drift on the water level through computer simulation results. To 
design the proposed controller the plant dynamics has been transfered into a suitable 
state space form. . 
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3.4.1 State Space Representation of dynamic model 

The dynamic equation (3.1) of pressurized water reactor for water level measurement 
can be splited into four state space equations as, 

xi(t) = Gi(g w (t) - q v (t)) 

x 2 (t) = -T 2 l x 2 {t) - (G 2 /r 2 )(q w (t) - q v (t)) 

x 3 (t) = -2Tf 1 x 3 (t) + Xi(t) + G 3 q w (t) 

X 4 ,(t) = -(T 1 -2 + 47r 2 T -2 )x 3 (f) 

Output water level is given by 

y p = xi(t) + x 2 (t) + x 3 (t) (3.2) 

The state space equation can be written in matrix form, 

Xp(t) = ApXp(tj T B p q w (t) + F p q v (t) 

y p (t) = C p x p (t) (3.3) 

Where, 

0 0 0 0 

0 — rT 1 0 0 

A> = 

0 0 -2 Tf 1 1 

0 0 — rf 2 + 47r 2 T _2 0 
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1 1 

Bp = Gi -G2T2 1 Gz 0 

1 T 

f p = I -G x G2T2 1 0 0 
c p= 1110 

-]. T 

x p(t) ~ Xi(t) X 2 (t) Xz{t) X t (t) 

The matrixes B v and F p can be written as 
B p = [b u 6 2j 6 3) 0] t 
F p = [dl,d2, 0, 0] T 

3.5 LQR Design for Steam Generator 

3.5.1 Introduction 

The linear quadratic regulator (LQR) is an optimal control problem where the state 
equation of the plant is linear, the cost function of the plant is quadratic and the 
test conditions consist of initial condition of the state and no disturbance inputs. The 
resulting optimal control known as the linear quadratic regulator. It can be applied 
to a wide range of control problems due to the flexibility available when selecting the 
quadratic cost functions [10]. The mathematical derivation of the LQR controller is 
given in appendix. 
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3.5.2 Problem Formulation 


For achieving asymptotically stable tracking with minimizing quadratic performance 
index by using LQR [11], the performance index is chosen as, 

Minimize J = J[Xp(t)Qx p (t) + Up(t)Ru p (t)]dt (3.4) 

subject to the state space equations as given in (3.3), 

x p (t) — A p x p ii^ 4” jB p Quj(t^ T FpQ'u(t) 

Vp{p') == C p x p (t) 

It is required that 

y p -+ y r as t -» oo (3.5) 

Here, y r = Water level reference. Its value depend on the design data of U tube steam 
generator. 

Tracking problem is converted into regulator problem by introducing a new state q 

q(t') = Vp Ur ~ C p x p y T (3-6) 

After introducing a new state in equations 3.3, the new state space equations in matrix 
form are given by 

X p q ApqXpq (t'j J3 p qUpq (t) T FpgU/(t) T H p qy r 

Vpqii) ~ CpqXpqii) (3.7) 
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where, 

control input u pq (t) = q W) 



- 


“ 

- 


- 


- 

Xpg 

Xp 

ii 

o 

Ap 

0 

, Bpq — 

B p 

F = 

> 1 PQ 

F P 


q 



o 

l 


1 

O 

i 


1 

o 

1 


a 


?? 


c p 0 


H pq = [0 0 0 0 - if. 


The cost function defined in equation (3.4) is thus reduced using equation (3.7), to 
the following form: 


Jl = J [Xpq(t)QlXpq(t) + v% q (t)Riu pq (t)]dt 

where Qi is state weight matrix and Ri is control weight matrix. To make R non 
negative, R\ should be positive definite and Qi should be positive definite or at least 
positive semidefinite. Ri and Qi are selected by the designer. For this particular plant, 
Ri and Qi are chosen to be unit matrices. 

The optimal control input is given by U pq (t) = —KX pq (t) 
where K is a constant feedback gain matrix given by 


K = R^B^Pq (3.8) 

where Pi is a constant symmetric positive definite matrix which is the solution of 
the Algebraic Riccati Equation (ARE) . The (ARE) is given below. 


P v-A-pq + Apq T Pi + Qi + PiB pq R 1 1 B pq Pi — 0 (3.9) 
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Time (sec) Time (sec) 

(a) water level (m) (b) water flow rate (kg/sec) 

Figure 3.2: LQ tracking at 50% Power 

The existence and uniqueness of solution for the above equation are guaranteed by 
the following assumptions: 

1 . (Apq,B pg ) is a controllable pair. 

2. A pq , Q 1 ^ 2 is an observable pair. 

3.5.3 Simulation Results 

A linear parameter varying model of UTSG of which parameters depends on the reactor 
power level. The model is linearized over four regions. Referring to table 3.1, the regions 
are divided according to operating power as: Region I for 0% < power < 15%, Region 
II for 15% < power < 30%, Region III for 30% < power < 50%, and Region IV for 50% 
< power < 100%. Over each region the elements of model matrices are assumed to 
vary linearly. The simulation studies are carried out for instances when the actual plant 
parameters drift from the linearized model parameters. It can be seen that variation of 
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Figure 3 . 5 : LQR Controller at 5 % Power at 3 % Perturbation 

element 62 of B p matrix has the most destabilizing effect on the system response. This is 
because it contains the term G2 and r 2 whose value determine the non-minimum phase 
characteristic of the steam generator. The simulation results are displayed in order 
to show the effects of these two parameters particularly on the effectiveness of the 
controller and the quality of the system response. As the power of operation becomes 
higher the effect of perturbation of parameters becomes less. It is assumed that water 
flow rate and steam flow rate are equal before the change in power demand occurs, that 
means the change in water level is zero. The simulation result are taken for the the 
power level of 5 % and 50 %. 

The simulation results are taken for the system response at the same power of 
operation ( 5 %) but the parameters (?2 and 72 are perturbed such that the element 62 
is increased by 3 % for 5 % power level and 30 % for 50 % power level. At steady state 
the water flow is matched with the steam flow which is desired. As shown from figures 
3.2, 3.3, 3.4 and 3.5 a substantial change in system behavior is noticed due to the drift 
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of system parameter, such as larger overshoot, more settling time and more oscillations 
before settling down to the desired value. Even when the parameter drifts, the controller 
is capable of maintaining the change in water level to zero and hence the water flow rate 
change equals the steam flow rate change. If the element b 2 drifts further the controller 
may not be robust enough to maintain the water level to the desired value. So we can 
say that the controller is robust enough to maintain the desired system response if the 
parameter perturbation lies within certain bounds. 

For comparison, we have displayed the simulation results. At 50% power the con- 
troller can tolerate maximum 30% perturbation. At 5% Power it can tolerate up to 
only 3% perturbation. 

3.6 FLC with GA Optimization 

3.6.1 Introduction 

In conventional control, the amount of control is determined in relation to a number of 
inputs using a set of equations to express the entire control process. Expressing human 
experience in the form of a mathematical formula is a very difficult task. Fuzzy logic 
provides a simple tool to interpret this experience into reality. In recent years, fuzzy 
logic has emerged as a powerful tool and is starting to be used in various applications. 
The application of fuzzy logic control technique appears to be most suitable one when- 
ever a well-d efin ed control objective cannot be specified and the system to be controlled 
is a complex one or its exact mathematical model is not available. Fuzzy Logic Con- 
trollers (FLCs) are robust and have relatively low computational requirements. They 
could be constructed easily using a simple microcomputer. 
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In present work a fuzzy logic controller optimized with genetic algorithm has been 
designed for controlling the water level of U tube steam generator. Input parameter for 
the fuzzification plant are 

e C0 = Vr - y P (t) (3.10) 

Now taking the derivative of the equation 3.10 we get 

e(t) = 0-y p (t) (3.11) 

where e(t) is the error in water level which is a function of time, 

Y r is reference water level in meter, which is a constant value. The value of y r depends 
on load. 

y p (t) is actual water level in meters. 

3.6.2 Fuzzification 

According to the fuzzy control strategy crisp inputs e(t) and e(t ) are converted to fuzzy 
inputs. For fuzzification, triangular membership functions are chosen with 9 regions. 
There are two inputs to the fuzzification plant. Total 81 rules are generated [12]. 
Maximum numbers of rules are fired for a given value of e, e are 4. 

3.6.3 Inference mechanism 

For generating a rule, two fuzzy sets are associated. Individual rule based inference 
mechanism has been used. In this mechanism the idea is to compute the overall value 
of the output variable based on the individual contributions of each rule in the rule 
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base. Each such individual contribution represents the value of the output variable as 
computed by the single rule only. In this problem, minimum of (/z(e)(t), pi(.e)(t) has 
been chosen as the truth value of a n th rule to be fired. 


3.6.4 Defuzzification Module 

It is described in section 2.7.3. In this problem, weighted average method [13] has been 
used for defuzzification. According to this, the control action is computed as, 


E n 

- i=l M 

Hw vr-\n 

2_/t= 1 Mi 


(3.12) 


where, 

g w =water flow rate ( control input). 
n=number of rules. 

/^membership for i t h rule. 
u — value of that particular rule. 

Finally this value has to be renormalized so that the point wise value of the control 
output is converted to its physical value. 

Genetic Algorithm for multi objective optimization (Nsga2 code) has been used 
for optimizing tool for the FL [14]. The parameters have been optimized are, 

1. Membership function for e(t) (error in water level) 

2. Membership function for e(t) (change in error in water level) 

3. Rules for q w (water flow rate) 

The following cost functions are opted to minimize the control input and error, 
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(a) Output (m) (b) Input (kg/sec) 


Figure 3.6: Water Level with FL optimized by GA at 50% Power 


t 


minimize'^ e 2 (t) 

i = 0 

t 

(3.13) 

minimize ^ 

z= 0 

(3.14) 


3.6.5 Simulation Results 

As described in section 3.5.3, the values of parameters are chosen accordingly. The FLC 
has been tested successfully at 50% of power and 50% of power with 3% perturbation 
of b 2 element of B p matrix . At steady state the water flow is matched with the steam 
flow which is desired. As shown from figures 3.6 and 3.7 a substantial change in system 
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Figure 3.7: Water Level of U Tube Steam Generator by FL optimized by GA at 50% 
Power at 3%perturbation 

behavior is noticed due to the drift of system parameters, such as more settling time 
before settling down to the desired value. Even when the parameter drifts, the fuzzy 
logic controller is capable of maintaining the change in water level to zero.. If the element 
&2 drifts further the controller may not be robust enough to maintain the water level to 
the desired value. So we can say that the controller is robust enough to maintain the 
desired system response if the parameter perturbation lies within certain bounds. 


Level of PWR 

Settling time in second 

Fuzzy controller with 50%power 

55 

Fuzzy controller with 50%power,3%perturbation 

75 

LQR controller with 50%power 

38 

LQR controller with 50%power,30%perturbation 

75 

LQR controller with 5%power 

300 

LQR controller with 5%power,3%perturbation 

320 


Table 3.2: Controller for Fuzzy and LQR Controller 
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3.7 Summary 


Control of UTSG water level strongly affects nuclear power plant availability. To at- 
tain a robust and stable control system, two controllers have been presented. LQR 
and fuzzy logic controllers have been implemented for UTSG water level control of a 
nuclear reactor. These controllers have been tested from 0% to 100% of the load. Both 
the controllers perform well when large-scale perturbation occurs. It is clear from table 
3.2 that the performance of LQR is better than the fuzzy logic controller in terms of 
settling time. It is shown that the tracking error can be reduced to zero if the param- 
eter perturbation is bounded within a certain value. Simulation results are provided 
for some cases to validate the effectiveness of the controller. After implementing the 
controllers, the water flow matches the steam flow, under various operating conditions 
and parameter uncertainties. The oscillations of water level and the settling time at 
50% power level are less than that at 5% power level, which implies that the control 
at low power is difficult. At high power operation, the performance of the controller 
is better. The controller enables the level to be kept within limits and the transient 
phenomena are overcome within reasonably short time. 



Chapter 4 


Design of a Steam Generator Level 
Control Loop for Pressurized Heavy 
Water Reactor 

4.1 Introduction 

Nuclear power generation is considered to be one of the major electric power source in 
future. The schematic diagram of nuclear power station is shown in figure 4.1. These 
power plants are very sensitive to any disturbances, (i.e, mechanical and electrical) 
due to various technical and safety reason. U tube steam generator is a major part of 
nuclear power production. It form a vital link between heat source (the reactor) and 
the heat sink (the turbo generator system) in most of the nuclear power plants. The 
schematic diagram of a U tube steam generator for of 500MW is shown in figure 4.2. 

Steam generator level control (S.G.L.C.) is designed to provide continuous heat sink 
to the reactor and ensures the reactor safety during the power production [15]. During 
reactor shutdown, the system requires removal of the heat. At any power level and 
under steady state condition, mass balance is achieved as well as level is maintained 
constant at desired set point. Due to any reason if the level control is lost the mass 
balance is affected. Level may rise and fall. If water level rises very high, moisture 
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Figure 4.1: U tube steam generator and its environment 

carry over into the turbine may occur. Hence turbine trip is proposed at very high 
level. If water level goes very low, it will affect heat removal and may affect reactor 
safety due to reduction in heat sink capability. For avoiding such situation, reactor step 
back is proposed at low level and reactor trip is proposed at very low level in the Steam 
generator. 

, In the process control loop, controller provides adjustable parameters which are 
adjusted to get the required system performance. If all the parameters are not set 
properly it may lead to poor stability and in the worst case instability. Study of the 
control loop by varying these parameters in a physical plant is usually not suitable due 
to various constraints. In simulation the physical system is replaced by its mathematical 
model. Response of the system towards various control schemes have been studied. 

The objective is to develop a simulation package for the steam generator level control 
loop for 500 MW PHWRs. It will help to study the transient behavior of the lopp and 
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aid in establishing the parameter values for optimum performance based on various 
performance indices. 

4.2 Description of the System 

Nuclear steam generator (SG) is a complex system. As shown in figure 4.2, the nuclear 
steam generator of 500MW PHWR is U tube type and consists of an integral tube 
bundle in a shell, moisture separator (MS), dryer and down comer region. Feed water 
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is pumped into the down comer region of SG through boiler feed pumps. Feed water is 
mixed with recirculated saturated water (coming from MS) and enter into heat exchange 
section. Heat rises in the shell around the tube carrying heavy water. Heat transfer 
takes place and the feed water gets heated up. Steam and water mixture is formed. 
This water goes to the moisture separator. Steam and water get separated and bulk 
of water is recirculated through the down comer. Steam along with remaining water 
(carry over) goes to dryer which further improves the steam quality. This steam is 
allowed to expand in the turbine or dumped directly to the condenser or vented to the 
atmosphere [16]. 

A simplified model of the steam generator water level is simulated. The following 
assumptions are. made for the analysis. 

(a) Conservation of mass. 

(b) Conservation of energy (enthalpy) . 

(c) Being a base load plant the heat input from the PHT system is constant. 

(d) Feed water temperature is constant. 

(e) Pressure in the system is an independent parameter and is uniform throughout 
the steam generator. Pressure is different at different electrical power. The profile 
of pressure with change in power is assumed to be known. 

The steam generator is divided into three parts: 

(a) Non boiling zone with a fixed boiling boundary. Temperature in this zone varies 
uniformly. 
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(b) From boiling boundary to average height of U tubes, steam quality varies linearly 
with height. 

(c) Above average height of U tubes to the level in MS (moisture separator), steam 
quality is assumed to be uniform. 


4.3 Factors Affecting the Level Change 

Change in level inside steam generator may occur because of change in 

(a) Feed water flow/Steam Flow. 

(b) Steam Generator Pressure. 

4.3.1 Effect of Feed Water Flow on Level Change 

At steady state, steam outflow and feed water inflow are same. Due to change in 
feed water inflow, there will be disparity between two flow rate and hence change in 
inventory over time. This will result in the level change. Since the enthalpy of feed 
water is less than that of saturated water. With the change in feed water flow, steam 
quality also changes. So it will change the homogeneous density and hence there will 
be change in water level. 

4.3.2 Effect of Pressure on Level Change 

Saturation temperature is a function of pressure. Hence the part of steam will be 
converted to water and vice versa, depending on the change in pressure from steady 
state condition. This will change the steam quality, homogeneous density and hence 
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level due to void collapse. Density of water, saturated steam and saturated water are 
function of pressure. So change in pressure will change the level directly. 

4.4 Model Development 

Based on the assumption of mass and enthalpy balance, a simplified model of steam 
generator level control is simulated in the following way. 

The change in steam quality due to change in feed water flow is calculated. The 
capacity of the system is accounted. Change in steam quality due to pressure change 
inside the steam generator is calculated. Overall effect of the above two changes are 
accumulated to find out the resultant steam quality. 

Due to change in steam quality, the volume change in the moisture separator is 
calculated. Finally the new level for steam generator is calculated considering volume 
changes. 

4.4.1 Different Parameters of this Model 

The following parameter values are taken for designing the steam generator of UTSG 
of 500MW Nuclear power plant. 

Internal shell diameter = 2.197m (meters) 

Internal Shell diameter in MS region = 3.35m 
External diameter of U tubes = 0.019m 
Number of U tubes = 2475 
U tube height in steam generator =10. 86m 
Steam generator moisture separator level = 13.09m 
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Internal diameter of moisture separator = 0.48m 
Number of moisture separators = 23 


-14.25m 

“14m 

“13.4m 

— 12.6m 


Turbine trip 


SG Level at 100% Load 
SG Level at 0-20% Loac 


Range of 
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of feed water 
control valve. 


Low 

level - . 

Alarm 


10.5m 


Reactor step back. 


V.lowl 
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Alarm! 


"7.85m 


Reactor trip. 


Figure 4.3: Control and Protection set point of Steam Generator Level 


4.4.2 Design specification 

Maximum pressure of the boiler=75kg/cra 2 
Pressure at 0% to 25% power =60kg/cm 2 
Pressure varies linearly from 25% to 100% of power 
Pressure at 100% power=41.9kg/cm 2 
Height of water Level for 0% to 20% of power=12.6m 
Height of water Level for 100% of power=13.4m 
Height varies linearly from 20% to 100% of power. 

The vaxiation of pressure water level of steam generator is shown in figure 4.4. 
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4.4.3 Transfer Function of the Process 

The process block is a part of the overall level control block as shown in figure 4.6. Water 
input is given to the steam generator and it has certain water level to be measured. The 
process of the steam generation is given in figure 4.5. As shown in this block diagram 
the whole process is divided into different blocks, 
steam flow rate. 

The block K\ is a constant gain which is converting the mismatch between water flow 
rate and steam flow rate to change in volume. 

The block K 2 is a constant gain which converts change in volume to to change in water 
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Figure 4.5: Process block diagram of steam generator 


level. 

The block K 3 is a constant gain which converts steam flow rate to the outlet electrical 
power. 

The different parameters used in this block diagram are, SV 5V \ and SV 2 the changes 
in volume (meter 3 ). 

6P is change in pressure 

S q , S gw , Sg P are the steam qualities. 

There are some terms related to process of steam generation as given in appendix. 
In the analysis of this process for the nuclear power plant the following two parameters 
are chosen to be constant. 

Temperature of feed water to the steam generator (Tp w )= 180 0 C. 

The dynamics of Block \ , Blocks and Blocks are described as below. 
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Mathematical Model of Block 1 


Inside the stem generator when the unsaturated steam is passed through the moisture 
separator the water is separated and recirculated . The temperature of the recirculated 
water is the saturation temperature of water. The recirculated water is mixed with feed 
water and converted to steam . 

Let the different parameters of the steam generator are written as, 

F ms — Water and steam flow rate in moisture separator (kg/sec). 

Sf = Steam flow rate (kg/sec). F w (t )= Feed water flow rate (kg/sec). 

Ryj(t )= Recirculated water flow rate (kg/sec) . 

W/(f)=Net water flow rate (kg/sec). 

T(t)=Temperature of net water flow rate (°C). 

7>=Temperature of feed water (°C ). 

Tr= T emperature of recirculated water flow (°C). It is saturation temperature of water. 
T s (t ) =Saturation temperature of feed water (°C'e). 

L(f)=Latent heat of vaporization of water (Joule/kg) 

E w (t)= Available energy to convert water to saturated water (Joule/sec). 

E s (t )= Available energy to convert water to unsaturated steam. 

13(f) = Available energy to convert saturated water to unsaturated steam (Joule/sec). 
S q (t )= Steam quality. 

5 9U) (t)=Steam quality due to change in water flow rate. 

5 , =Specific heat of water. It is taken as as 4.18 kilo Joule/kg in this problem. 
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The value of the above mentioned parameters are calculated as 


q ^ Water + Steam 
Sf(t) = F ms X Sq(t) 

Rw(t) = F ms X (1 - S q (t)) 

Wj(t) = Rv(t) + F w (t) 


The temperature of water after recirculation is given by, 


Tt(t) = 


F w (t) x Tp + Rw(t) x Tr 
F w (t ) + Ftw(t) 


(4.1) 


Heat exchange takes place from primary side to secondary side of U tube steam gener- 
ator. It is utilized to convert water into steam. In each iteration there is a change in 
water flow rate. So there is a change in energy of the water. Finally it results a change 
in steam quality [17]. 

The energy gained by water for converting water to unsaturated steam is given by, 


E s (t - 1) = W f (t - 1) x (S x ( T,(t - 1) - Ti(t - 1)) + L{t - 1) x S q {t - 1)) (4.2) 


The energy gained by water for converting water to saturated water is given by 

E w (t) = W f (t) xSx (T s (t) - Ti(t)) (4.3) 

The energy available for converting water into steam can be calculated from equations 
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(4.2) and (4.3) as, 


E(t) — E a (t — 1) — E w (t) 


(4.4) 


Energy E(t) is not sufficient for converting all the water inside the steam generator into 
steam. Prom equation (4.4), steam quality due to water flow rate can be written as 


S gw (t) = 


E 

W f (t) x L[t) 


(4.5) 


Prom equation (4.5) it can be conclude that the steam quality changes due to change 
in water flow. 


Mathematical Model of Block 2 

Due to change in pressure, there is a change in energy in the boiler. Change in energy 
is responsible for change in steam quality. So there is a change in steam quality due to 
change in pressure. 

As described in earlier section, the steam generator is divided into different zones. 
The mass in the different zones are, 

M) >u = Total mass inside the steam generator from boiling section to U tube section 

(kg)- 

In this zone both saturated water and saturated steam are present. Steam quality varies 
linearly from zero to maximum value ( S g ). The average value of steam quality is chosen 
for simulation. 

Mu m s = Total mass from U-tube to moisture separator (kg) from U-tube to moisture 
separator. In this zone Steam quality is assumed to be uniform. 

M ma = Total mass from moisture separator to steam outlet (kg), in this zone steam 
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quality is assumed to be uniform . 

AE p (f)=Extra energy of water due to change in pressure. (Joule/sec) 

E p (t)~ Energy of the steam (Joule/sec). 

E ps (t)= Energy available for converting saturated water to unsaturated steam (Joule /sec). 
Sqp{t )= New steam quality due to change in pressure. 

The change in energy of water inside the steam generator due to change in pressure 
is given by, 


E p (t—1) = (Mb u (t—l)x(l—S q (t—l)/2)+(Mu rns (t—l)+M ms (t))x(l—S q ))x(SxT sp (t—l)) 

(4.6a) 

E p (t) = ( M bu (t ) x (l-S q (t)/2) + (M Um ,(t) + M ms (t)) x (1-5,)) x (SxT sp (t)) (4.6b) 

A E p (t) = E p (t - 1) - E p (t). (4.6c) 


The energy of the steam can be written as, 


E p (t) = (M bu (t)/2 + M Uma {t) + M ms (t )) x L(t) x S q (t) (4.7) 


From equations 4.6 and 4.7, the energy for converting saturated water into steam is 
derived as. 


Ep S (t ) = A E p (t) + Ep{t ) 


(4.8) 
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Now, the steam quality is found to be 


S (t) = E ps (t) 

qA } (M bu (t)/2+ M Um ,(t) + M ms (t)) x L(t) (4 ' yj 

Mathematical Model of Block 3 

According to the description given in section 4.4.3, there is a change in steam quality 
due to change in pressure. This results in change the homogeneous density in the boiler 
as well as change in volume of water. 

The different parameters for calculating the process are given below, 
D sw (t)=Density of saturated water (N/m 2 ). 

A, s (t)=Density of saturated steam (N/m 2 ). 

Duh(t)= Homogeneous density from boiling section to U tube section (N/m 2 ). 
DMSh{t)— Homogeneous density from U tube section to steam outlet (N/m 2 ). 

Al/i= Volume of water reduced from boiling section to U tube section in meter 3 . 
A^=Area of cross section in U tube section in meter 2 

Hu — H b b= Height of steam generator from boiling zone to U tube zone (meter). 

A V 2 = Volume of water reduced from U tube to moisture separator region in meter 3 . 
A 0[/ =Area of cross section in moisture separator (meter 2 ). 

Hms — Hu = Height of steam generator in moisture separator in meter. 

AV 3 =Volume of water reduced from moisture to steam outlet region in meter 3 . 
A ms =Area of cross section above moisture separator to steam outlet in meter 2 . 
lev= Height of steam generator above moisture separator. 

AV=Total volume reduced due to change in steam quality. 

The homogeneous density of steam and water from boiling section to U tube section 



due to pressure in previous iteration is given by, 


Duh{t) 


Duh{t — 1) 


S qP (t)/( 2.0 x £>„(*)) + (1 - S qp (t)/2.0)/D sw (t) 

1 

(S qp (t - 1)/ (2.0 x D u (t ~ 1)) + (1 - S gp (t - l)/2.0)/D sw (t - 1))' 

DMSh.it ) — 77; TTwn TTTTTi cv TTvTTri 777 ( 4 - 12 ) 


(4.10) 

(4.11) 


DMSh{t — 1) = 


C Sq P (t)/D ss (t ) + (1 - St qp (t))/D sw (t) 

1 

(Sji - 1)/D sa (t - 1) + (1 - S^t - l))/£> w (i - 1)) 


(4.13) 


Decrease in volume due to change in homogeneous density, 


A Vi = A v x ( H v - H bb ) x 
AV 2 = A au x ( Hms ~ H v ) x 
AV 3 = A ms x Lev x 


D_ uhif) — Duhjt — 1 ) 

Du hit) 

Dushjt) - PMShi t — 1 ) 

DMSh(t ) 

(Dm shit) — DMSh(t — 1 )) 


Dmsh (t) 

AV = A Vi + AV 2 + AV 3 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


Prom equation (4.17) the volume change due to change in pressure is calculated. 

Block of look of table 

This block convert electrical power to to pressure according to design specification as 
shown in figure 4.4 


4.5 Steam Generator Level Control Scheme (SGLCS) 

For steam generator (SG) level control, three elements are controlled, i.e., water level, 
feed water flow rate and steam outflow from the S.G., 
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With the increase in load, steam flow will increase . So there will be more steam 
bubble in SG, which will cause water to swell or raise rather than fall. Thus there will 
be a transient increase in the water level. Similarly when steam flow reduces or feed 
water flow increases, bubble tend to collapse and temporarily reduce their formation. 
This will cause water level to fall rather than rise because of reduced water usage or 
increased water supply. This will drive the control signal in wrong direction, if control 
is done with level signal alone. This phenomenon is called “Inverse Response” . It will 
cause delay in control action which makes the controller design more difficult . 

Steam generator level is programmed to increase linearly from no load to full load. 
The Steam generator level program is shown in figure 4.3. The set point for water level 
in steam generator has been fixed for safety and protection of nuclear power plant. It 
varies linearly from 20% of load to full load. The controller has to be designed in such 
a way that the water level in steam generator should follow the set point in order to 
avoid the reactor shut down or back stepping of the reactor. 

The level is measured by the level transmitters as shown in figure 4.6. Its value 
is compared with fixed level set point. Resulting error is fed to the controller. From 
the controller output, signal from the feed water flow is subtracted and fed to the 
control valve through E/A (electrical to analog) converter. Control valve is modulated 
to control the flow of water under the steady state and hence the water level in SG is 
controlled to the designed value. 

As the transfer function of the process, i.e., steam generator is very much compli- 
cated, the formulation requires special attention. As is shown is shown in the figure 
4.6 the transfer functions of various other elements in the control loop, i.e, E/A, con- 
trol valve (CV) converter, flow transmitter, level transmitter etc are very simple. The 
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coefficient of the transfer function is worked out for the design information according 
to the manufacture data. The overall block diagram of SGWL is shown in figure 4.6. 
The various blocks are associated with it are “ primary controller, process, secondary 
controller, E/A (electrical to analog converter), CV (Control valve), FLOW Tx, Level 
Tx”. The transfer function of each block diagram are described below. 


X l=Freed water flow rate 



Secondary Controller 

As given in the block diagram of steam generator level control in figure 4.6, The sec- 
ondary controller is chosen as a proportional controller. So its transfer function is, 


TF = K p i 


(4.18) 
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Where 


K p i=Gain of the Controller. 

E/A Converter (E/A) 

Basically electrical to analog (E/A) converter is a delay circuit. There is some time 
delay in converting the electrical signal to analog signal. So the transfer function is, 


TF = 


G 

1 + sT 


(4.19) 


where, 

Gain (G) =1, time constant (T)=0.2 second. 


Control Valve (CV) 


The control valve is used to open the valve of feed water flow. The transfer function 
for control valve is 


TF = 


G 

1 + sT 


(4.20) 


where, 

Gain (G)=l, time constant (T)=0.4 second. 


Feed Water Flow Transmitter (FLOW Tx) 

It is basically a transducer which convert feed water flow to the equivalent electrical 
signal. The output of the flow transmitter is fed to the secondary controller. The 


55 



transfer function of the feed water flow transmitter is 


TF = 


G 

1 + sT 


(4.21) 


where, 

Gain (G)=l, time constant (T)=0.22 second. 

Water Level Transmitter (Level Tx) 

It is basically a transducer which converts water level to the equivalent electrical signal. 
The output of the Tx can be used as as a feed back signal to the primary controller. 
The transfer function of the feed water flow transmitter is 

TF =iT7f (422) 


where, 

Gain (G)=l, time constant (T)=0.2 second. 

4.5.1 Controller Design 

As shown in the figure 4.6, a cascade control scheme has been used for water level 
control of nuclear steam generator. It consists of two controllers. The inner control 
loop is called secondary control loop and the outer control loop is called primary control 
loop. The secondary controller is chosen as a proportional controller that is derived for 
the given problem. The primary controller is chosen as proportional derivative (PD) , 
proportional integral(PI) and fuzzy logic (FLC) controller. 


56 



4.6 Proportional Derivative Controller 


A proportional derivative controller has been implemented as a primary controller for 
the given problem. The error in water level to the primary control loop is given by, 

e(t) = y r ~ Va{t) (4.23) 


where , 

y r =Reference water level in meter, 
yjt) is the actual water level in meter. 
Taking its time derivative, we get 


e(t) = 0 - y a (t) 

The control output from the proportional derivative controller is 

control = K p x e(f) + K d x e(t) 


(4.24) 


(4.25) 


where, 

K p = Proportional gain of primary controller. 

Kd= Derivative gain of primary controller. 

K p \= Proportional gain of secondary controller. 

The optimal values of K p , Kd and K p i are optimized by using genetic algorithm. 
Following cost functions are chosen to minimize the “water flow rate” and “error in 
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water level” . 


t 


minimize e 2 (f) 

i=0 

+ 

(4.26) 

minimize ^ q^t) 

i=0 

(4.27) 


where, 

q w is the water flow rate in kg/sec. 

After optimization the value of K p , Kd and K p \ are found to be 599.299622, 80.664223 
and 1.191765 respectively. 
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(b) Inlet Water Flow 


Figure 4.7: Water Level Controller by Using PD Controller 


4.6.1 Simulation Results 

The proportional derivative controller optimized with genetic algorithm is implemented 
for controlling the water level of U tube steam generator. The simulation results are 
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Figure 4.8: Water flow and steam flow rate for PD control 

in the figure 4.7. It is observed that, there is no overshoot but some steady state 
error remains in the desired output. After optimizing the values of parameters for 
proportional derivative controller gain, there is an improvement of the performance, but 
the steady state error can not be eliminated completely. Under steady state condition 
the water flow is matched with the steam flow which is desired. It has been shown in 
figure 4.8 


4.7 Proportional Integral (PI) Controller 

A proportional integral (PI) controller has been implemented as a primary controller 
for the given problem. As described in section 4.6 the error in water level to the primary 
control loop is given by, 

e(t) =y r - y«(t) (4.28) 
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The control output from the PI controller is, 


controli = K p x e(f) + f Ki x e(f) (4.29) 

JO 


where, 

K p = Proportional gain. 

Ki= Integral gain. 

The optimal values of K p , Ki and K p i are calculated by using genetic algorithm. 
Following cost functions are chosen to minimize the “water flow rate” and “error in 
water level” during the simulation. 


t 


minimize e 2 (t) 

i=0 

+ 

(4.30) 

minimize g^(f) 

i=0 

(4.31) 


After optimization, the values of K p , Ki and K p i are chosen to be 249.413925 ,49.970673 
and 1.162353 respectively. This output of PI controller is fed as the reference to the 
secondary control loop. 


4.7.1 Simulation Results 

The PI controller optimized with genetic algorithm is implemented for controlling the 
water level of U tube steam generator. The simulation results are shown in the figure 
4.9. The set point tracking is achieved by using this controller. It is observed that, 
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(a) Water Level (b) Inlet Water Flow 

Figure 4.9: Water Level Controller by Using PI Controller 



Figure 4.10: Water flow and steam flow rate for PI control 

initially there is some overshoot but it damped out after some time. Under steady state 
condition the water flow is matched with the steam flow which is desired. It is shown 
in figure 4.10 
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4.8 Fuzzy Logic Controller 


A fuzzy logic controller (FLC) is implemented as a primary controller for the given 
problem. The mechanism of the fuzzy logic controller is described in chapter 2. As 
described in section 4.6 the error in water level to the primary control loop is given by, 

e(i) =y r - y a (t ) (4.32) 

As described in earlier section, fuzzy logic controller is applied to this plant. The pa- 
rameters of fuzzy logic controller such as “membership function for error” , “membership 
function for change in error”, “rules” are optimized by using genetic algorithm. Fol- 
lowing cost functions are opted to minimize the “water flow rate” and “error in water 
level” through out the simulation. 


t 

minimize yym (4.33) 

t= 0 
t 

minimize X>’w (4.34) 

1=0 


After fuzzification and defuzzification the control output is generated from the fuzzy 
plant, is fed as a input to the secondary controller. 


4.8.1 Simulation Results 


The fuzzy controller optimized with genetic algorithm has been implemented success- 
fully for controlling the water level of U tube steam generator. The simulation results 
are shown in the figure 4.11. It is observed that, there is no overshoot and also steady 
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Figure 4.11: Water Level Controller by Using Fuzzy Controller 



Time (Seconds) 


Figure 4.12: Water flow and steam flow rate for FL control 

state error is approximately zero . The settling time is least among all the three con- 
trollers. Under steady state condition the water flow is matched with the steam flow 
which is desired. It is shown in figure 4.12 
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Level of UTSG 

Steady state error in % 

Overshoot in % 

Settling time in sec 

Fuzzy controller 


0 

1.5 

PI controller 


11.11 

20 

PD controller 

2.41 

0 

2 


Table 4.1: Settling Time for Different Controllers 


Description of the table 

The simulation results of three different controllers have been shown in table 4.1. It 
shows that for FLC controller, steady state error in % and settling time in second 
for water level are 0.06 and 1.5 respectively. For PI controller, steady state error in 
%, overshoot in % and settling time in second for water level are 0.03, 11.11 and 20 
respectively. For PD controller, steady state error in % and settling time in second for 
water level are 2.41 and 2 respectively. It is seen that the performance of fuzzy logic 
controller is better than other two controllers. 

4.9 Validity of the Model 

The described scheme is used for developing a nonlinear model of steam generator 
level control loop for 500 MW PHWR . But the geometry of the steam generator 
is also different from the one used for 500MW PHWR. Besides this there are some 
minor differences in designing the steam generator which are ignored in developing the 
steam generator level model. Parameters are obtained for theoretical observation under 
optimum performance where as the parameters obtained from practical observations are 
operating parameters which may not be optimum. There may be some errors in the 
designed model due to some wrong approximations. 

So, the system may try to operate around the theoretical optimum parameters first 


64 






by reducing the gain and then by slowly reducing the reset time to further validate the 
simulation. 

4.10 Summary 

A nonlinear model of steam generator level control loop for 500 MW PHWR has been 
proposed. The two controller are used in cascade mode so that full control capacity 
of the control valve can be utilized. In this scheme the primary controller is placed 
in the outer loop to control the water level in SG. Output of this controller is set 
point for secondary controller (Feed- water flow controller). Secondary controller is a 
proportional controller. The primary controller has been implemented by using fuzzy 
controller, PI controller and PD controller. By using above controllers, the water level 
of nuclear steam generator is controlled. The results of all the controllers have been 
shown in Table 4.1. It is observed that the settling time in case of fuzzy logic controller 
is smallest than other two controller. In PI controller, steady state error is zero but 
there is over shoot. In PD controller though there is zero overshoot but some steady 
state error is remained. But in Fuzzy logic controller there is no over shoot and steady 
state error is approximates to zero. 
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Chapter 5 

Nonlinear Control System Design: 
Single Link Manipulator 

5.1 Introduction 

In recent years considerable amount of interest has been paid by the researchers on 
the use of fuzzy logic-based controllers. Many successful applications have done and 
fuzzy logic controller are commercially available. Processes which are it is already 
commercially available. Processes which are controllable by a skilled operator by using, 
fuzzy logic controller rather than the use of conventional control methods are the target. 
Single link manipulator is a nonlinear single input single output system. It has one 
actuator and one degree of freedom. A flexible arm is connected to a motor. The 
angular motion of the arm depends upon the angular motion of the motor. The input 
torque is given to the motor which is under our control. So by controlling the torque, 
the position of the arm can be controlled. 

Industrial manipulators are often equipped with conventional PID controllers due 
to their simplicity in structure and ease of design. However, when using a PID control, 
it is difficult to achieve a desired level of control performance. One of the most widely 
used design methods for FLC is to define membership functions of linguistic variables 
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and to formulate fuzzy rules by control engineers [18]. Design of an efficient fuzzy logic 
controller involves the optimization of parameters of fuzzy sets and proper choice of 
the rule base. The parameters of the controller have been optimized based on genetic 
algorithm. 

The objective of the present research is to study more advanced control techniques 

for a single link manipulator in order to achieve a better performance both in simulation 

| 

and experimentation. Specifically, a fuzzy logic based control configuration is developed 
for precise position tracking control of a single link manipulator, various conventional 
control strategies have been implemented for position tracking. 



Y 

Figure 5.1: Single Link Manipulator 
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5.2 Model Development 


The single link manipulator is shown in figure 5.1. A arm is connected to the motor. 
The motion of the arm depend upon the motion of the motor. A dynamic model of the 
manipulator has derived by using Lagrangian formulation of equations of motion. The 
parameters have been to describe its motion are, 
m=Mass of the pendulum in kg. 

£=Length of the pendulum in meter. 

8— Angle of the pendulum from vertical bottom position in rad. 
g=Acceleration due to gravity in . 
r= Control torque applied to the pendulum in Nm. 


5.2.1 Lagrangian formulation of equations of motion 

The dynamic equation of single link manipulator has been derived in this section [11]. 
Prom the equation of motion, the kinetic energy K of the connecting rod is, 


K = \ml 2 6 2 

Zj 


( 5 . 1 ) 


We assume that the potential energy of the pendulum is zero when it is at rest at 
vertical downward position. Hence the potential energy U of the pendulum at any 
position is given by the equation, 


U = mgl{ 1 — cosd) 


( 5 . 2 ) 
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The Lagrangian function L is 


II 

1 

(5.3) 

ml 2 9 2 — mgl( 1 — cos9) 

(5.4) 


Now the Lagrange equation describing the pendulum motion is 


d_ dL _ 9L 
dVdd 96 


(5.5) 


Using the equations (5.4) and (5.5), the dynamics of single link manipulator (SLM) can 
be written as, 

ml 2 9 = —mglsind + r (5.6) 


The parameters values are used for performing the simulation are given below 
m — 1 kg, 1=1 meter, g = 9.81kg /meter 2 . 


5.3 Control Objective 

The pendulum behaves in to and fro motion which produce sinusoidal oscillation. The 
control force is generated in such a way that the motion of the pendulum follows si- 
nusoidal trajectory. The block diagram of the control system is shown in Figure 5.2. 
A reference sine wave is compared with actual position of the single link manipulator 
(SLM) given as the error in the plant. Different controllers such as proportional deriva- 
tive, feed back linearization, fuzzy logic controller which have been used are described 
below. 
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5.4 Proportional Derivative Controller 


A proportional derivative controller has been designed for the given problem. The error 
signal is, 

e(t) = sin(t ) — 9{t ) (5.7) 

t is the the time of the response in seconds. 

e(t) is the input to the controller. Now taking its time derivative, we get 

e(f) = cos{t) — 9(t) (5.8) 

The control torque for proportional derivative controller is 

t = K p x e(t) + K d x e(t) (5.9) 


where, 

K p — Proportional gain. 

Kd= Derivative gain. 

The optimal value of parameters K p and K d are calculated by using genetic algorithm. 
For applying genetic algorithm to the given system the following cost functions are 
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chosen, 


t 

minimize E e T (t)kie(t ) + e T k2e(t ) 

i — 0 


t 

minimize E t2 « 
i=0 


(5.10) 

(5.11) 


k \ , k% are taken, as unity. After optimization, the values of K p and are found to be 
74 and 3 respectively. 



(a) Pendulum Angle (rad) (b) Input (Nm) 


Figure 5.3: PD Controller for Single Link Manipulator 


5.4.1 Simulation Results 

The proportional derivative controller optimized with genetic algorithm has been im- 
plemented for designing the controller for single link manipulator. The simulation for 
this system has been done in C++ language. The simulation results have been shown 
in the figure 5.3. It is observed that, the desired trajectory of SLM does not follow 


71 




perfectly the path of sinusoidal response. There is an steady state error in the output. 
Initially there is some oscillation but latter it has damped out. 

5.5 Feedback Linearization Technique 

Feedback linearization is nonlinear control system design procedure which has attracted 
many researchers in recent years. The central idea of this approach is to algebraically 
transfer a non linear system dynamics into a (fully or partly) linear one, so that the 
linear control techniques can be applied. It can be used as model simplifying de- 
vices in development of robust or adaptive nonlinear controllers [19]. It has been used 
successfully to address some practical control problems. These include the control of 
helicopters, high performance air crafts, industrial robots and biomedical devices. 

In this chapter feedback linearization technique is used to designing the controller 
for single link manipulator. The details of feed back linearization is given in appendix 
The dynamic equation of single link manipulator, given in equation (5.6). It can be 
written as 

ml 2 6 + mglsin(6 ) = r (5.12) 

All the symbols are defined in section 5.2. Choosing sin(t) as trajectory, the error can 
be written as 

e(t) = sin(t ) — 9(t) (5.13) 
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Taking time derivative of the above equation we get, 


e(f) = cos(t ) — 9(t) 

(5.14) 

e(t) = — sin(t ) — 9(t) 

(5.15) 


A function “r(t) n is chosen such that 


r(t) = e(t) + A e(t) (5.16) 

where, A is a constant multiplying factor. 

Taking derivative on both sides of equation (5.16), 

r(f) = e(t) + A e(t) (5-17) 

Putting the value of e(t) from equation (5.14) in equation (5.17) we get 


r(f ) = — sin(t ) — 9(t) + A e(t) 


(5.18) 


Putting the value of 9(t ) from'equation (5.12) in equation (5.18) we get 


r(t) = — sin(t ) + A e(t) 


(r — mglsin(6 )) 

ml 2 


(5.19) 


Taking the values of m=l kg, and 1=1 meter, the equation becomes 


r(t) = -sin(t) + \e(t) - (r - gsin(e)) (5.20) 
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r is proposed in such a way that 


rft) = —krft) 

(5.21) 

Where A; is a positive constant. The solution of equation (5.22) 

can be expressed in 

time domain as 


rft) = e~ kt 

(5.22) 

from equation (5.22) it can be said that, 


lim rft) = 0 

t— >oo 

(5.23) 

=» lim (eft) 4- A eft)) = 0 

t-+ OO 

(5.24) 


From equation (5.24) it can be concluded that 

lim eft) — 0 (5.25) 

t— ► 00 

5.5.1 Simulation Results 

The feedback linearization controller has been implemented for designing the controller 
for single link manipulator. The feed back gain matrix is calculated by using Matlab. 
The simulation results are shown in the Figure 5.4. It is observed that the trajectory 
of SLM perfectly follow the path of sinusoidal response. There is no oscillation in the 
desired trajectory, control action is smooth. The performance of the controller is given 
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Figure 5.4: Trajectory tracking of SLM by FBL Controller 


in table 5.1. 


5.6 Fuzzy Controller with GA Optimization 

In present work fuzzy logic controller has been designed for sinusoidal tracking of the 
single link manipulator. Input to the fuzzification plant are taken from error in position 
and velocity (derivative of position). So it also called as fuzzy PD controller. 


5.6.1 Fuzzification 

According to the fuzzy control strategy, crisp inputs e(t) and e(t) are converted to 
fuzzy inputs. For fuzzification, triangular membership functions are chosen with 9 
region. There are two inputs to the fuzzification module. Total 81 rules are generated 
[12]. Maximum number of rules are fired for a given value of e(t), e(t) are 4. 
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5.6.2 Inference mechanism 

As described in section 3.6.3, individual rule based inference mechanism is used. In this 
problem, following inference has been chosen. 

Minimum of (/j,(e)(t), fi(e)(t)) is taken as the truth value of n th rule to be fired. 


5.6.3 Defuzzification Module 

It is described in section 2.7.3. In this problem, weighted average method [13] has been 
used for defuzzification. According to this, the control action is computed as, 


_ _ EILi 

n.1* 


(5.26) 


where, 

r=Control input. 
n=number of rules. 

/i=membership for each rule. 
o>=value of that particular rule. 

Genetic Algorithm for multi objective optimization (Nsga2 code) has been used as 
the optimizing tool for FLC [14]. The parameters have been optimized are, 

1. Membership function for e[t) (error in position). 

2. Membership Function for e{t) (error in velocity). 

3. Number of rules. 
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The following cost functions are opted to minimize the torque and error throughout the 
simulation, 


t 

minimize E e T (t)k 1 e(t) + e T (t)k 2 e(t) 
2=0 


t 

minimize r 2 
i = o 


(0 


The values of /c 1; have been taken as 3, and 1 respectively. 


(5.27) 

(5.28) 



(a) Pendulum Angle (rad) 



Figure 5.5: FL of Single Link Manipulator 


5.6.4 Simulation Results 

The simulation results are shown in the figures 5.6 and 5.5. It is observed that the 
trajectory of SLM follows the path of sinusoidal response. Control action is smooth. 
Figure 5.5 shows the trajectory without optimization and figure 5.6 shows the trajec- 
tory with optimization. From figures 5.5 and 5.6, it cleared that after optimization 
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Figure 5.6: FL with GA Optimization for Single Link Manipulator 


of certain parameters of FLC, the controller gives a better performance. Compara- 
tive performance of all the control techniques is given in table 5.1. Table 5.1 shows 


Single Link Manipulator 

RMS error of position in °C 

RMS Control force in Nm 

Fuzzy controller 

0.006818 

5.16258 

PD controller 

0.064365 

4.8038 

Linearized feedback controller 

0.0001125 

5.1755 


Table 5.1: Trajectory Tracking of SLM 


a comparative performance between FLC, PD, and PI Linearized feedback controllers 
in terms of RMS value of error and control force. It observed from this table that for 
FLC controller, the root mean square of error of position of pendulum in radian and 
root mean square of control force in Nm are 0.006818 and 5.16258 respectively. For PD 
controller, the root mean square of error of position of pendulum in radian and root 
mean square of control force in Nm are 0.064365 and 4.8038 respectively. For linearized 
feedback controller, the root mean square of error of position of pendulum in radian 
and root mean square of control force in Nm are 0.064365 and 4.8038 respectively. 
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5.7 Summary 


The system dynamics of single link manipulator (SLM) as derived using Euler-Lagrange 
formulation. It has a nonlinear dynamic. Three controllers such as fuzzy logic, pro- 
portional derivative, feedback linearized controller have been implemented for getting 
the sinusoidal trajectory of single link manipulator. Fuzzy logic controller optimized 
with GA gives a very good trajectory tracking. The control action is less and smobth. 
This gives us the motivation for using fuzzy controller for more complex non linear 
systems. Fuzzy logic controller is relatively easy to develop and robust as compared to 
conventional control strategies. 
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Chapter 6 

Nonlinear Control System Design: 
Reaction Wheel Pendulum 

6.1 Introduction 

The Reaction Wheel Pendulum (RWP) is a non linear under actuated system. It is 
a pendulum with a rotating wheel at the end, which is free to rotate about an axis 
parallel to the axis of the pendulum. The wheel is actuated by a dc motor, while the 
pendulum is under actuated [20]. The motor produces a torque on the wheel, causing 
the wheel to spin. The control torque is generated by the angular acceleration of the 
disk (pendulum) which is responsible to control the system. 

Balancing the pendulum in the vertical upward position is difficult. By using fuzzy 
logic and state feed back controllers, balancing controllers have been designed in this 
chapter. 

6.2 Model of the Reaction wheel pendulum 

m i : mass of the pendulum in Kg. 
m 2 : mass of the wheel in Kg. 
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Figure 6.1: Reaction Wheel Pendulum 
h : length of the pendulum in m (meter). 

l c i : distance from the pivot to the center of the pendulum in m. 
I\ : moment of inertia of the pendulum in Kgm 2 . 
h : moment of inertia of the wheel in Kgm 2 . 
qx : pendulum angle with vertical axis in radian. 

<7i= Velocity of pendulum in radian/sec. 

<?2 : angle of the wheel in radian. 
i?2= Velocity of the wheel in radian/sec. 

r : motor torque input (control input) applied on the disk in Nm. 
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6.2.1 Lagrangian formulation of equations of motion 

The kinetic energy of the connecting rod alone, is given by 

Ki = ^(mi& + m 2 l\)qi (6.1) 

The kinetic energy of the wheel is given by 

K 2 = 2 m 2^i9i 2 + 2 ^ 2 ^' 1 + q 2 ) 2 ( 6 . 2 ) 

Prom equations (6.1) and (6.2), the kinetic energy of the system is K = K\ + K 2 . 
where, 

K = -(mi l 2 x + m 2 l\ + 1\ + I 2 )qi 2 + I 2 qiq 2 + ^I 2 q 2 2 (6.3) 

Total potential energy of the system comes out to be 

P = (miUg + m 2 lig)(cos(qi) - 1). (6.4) 

Lagrangian equation of the system is given by L = K - P. 

L = -(mi/^ 4- m 2 l\ + Ii 4- I 2 )qi 2 + I 2 qiq 2 + -I 2 q 2 2 — m(cos(qi) — 1) (6.5) 

where m = mil c i + m 2 li 
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Prom Euler-Lagrange’s equations of motion 


d.dL dL 

iM {q ' q)) -T q {q ’ q) = T 


( 6 . 6 ) 


where, 

(^position in radian, 
g=velocity radian/sec, 
r=control input torque in Nm. 

Prom the equations (6.5) and (6.6), the dynamic equation of reaction wheel pendulum 


(mil 2 cl + m 2 l{ + h + h)qi + hq-i - mgsin(qi) = 0 

I 2 qi + I 2 q 2 = t 


(6.7) 

(6.8) 


In compact form, the dynamic equation of the system can be written as 


D(q)q + g{q ) = u 


(6.9) 


where, 





" 

q = 

q i 

, u = 

0 


q 2 


r 


D(q) is the inertia matrix and is given by 


D{q) = 


m l^cl d" ttl 2 1\ + I\ + 1 2 I 2 

h h 


“1 

” 

“ 


dn 

di2 


i 

tO 

1 
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9(Q) = 


—mgsin(qi) 

0 

Hence the dynamic equation can be written as 


duq'i + di 2#2 - mgsin(qi) = 0 

(6110) 

d2lil + d 2 2?2 = T 

(6.11) 


Pendulum parameters: 

The kit for doing real time experiment is available in intelligent control lab. The values 
of various parameter of the kit have been taken for the simulation of the system. 

Mass of the pendulum, mi = 0.2164 Kg 
Mass of the wheel, m 2 = 0.0850 Kg 
Length of the pendulum, li = 0.1270 meter 

Distance from pivot to the center of mass of the pendulum, l c i = 0.1173 m 
Moment of inertia of the pendulum, I\ =2.233 xl0~ 4 Kgm 2 
Moment of inertia of the wheel, / 2 =2.495xl0 -5 Kgm 2 
Ratings of the drive motor: 

Maximum speed of the motor, N max — 825 rad/sec. 

Electrical Time constant, Te =0.5 ms (milli second). 
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6.3 State Feedback Controller by Pole Placement 

The non linear dyanamics (6.10) and (6.11) are linearized around the unstable equlib- 
rium point which is vertical upward position of the pendulum. A state feedback con- 
troller has been designed using this linearized model. 

The mathmatical derivation of state feedback controller is given in appendix. The 
system dynamic (6.10) and (6.11) can be written in the following form: 


9 i = 
92 = _ 


d^mg . 

stn(q 1 ) 


d 


12 


det(D) 

-d 2 img 

det(D) 


det(D) 


SinWl) + ^D) T 


( 6 . 12 ) 

(6.13) 


where det(D) = d n d 2 2 — d^d^i- 


The linearized state space equations for the given plant is given as 


X = AX + Bt 


(6.14) 


X 


<?1 ?1 92 92 


A = 


“ 





- 

0 

0 

1 

0 


0 

0 

d22 mg 
det(D) 

0 

0 

, B = 

12 

det(D) 

0 

0 

0 

1 


0 

0 

-chimg 

0 

0 


dn 

det(D) 


det(D) 
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Where X is the state vector, A is the system matrix, B is the conreol matrix and r is 
the control torque. 

The controllability matrix Q for this system is given by [ B , AB , A 2 B, A 3 B ]. The 
rank of the controllable matrix is found to be 4 [21]. Hence the linearized system is 
controllable. A full state feedback controller has been designed. 

The control law r = —KX is chosen so as to place the poles of A — BK in the left 
half of s — plane. Then the system becomes asymptotically stable, 
where K=[K U K 2 , K 3 , K 4 }. 




Figure 6.2: SFB Controller for Reaction Wheel Pendulum 
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Figure 6.3: Control Force for SFB Controller 

6.3.1 Simulation Results 

The state feedback (SFB) controller has been implemented for reaction wheel pendulum 
which is a nonlinear system. The simulation have been shown in the figures 6.2 and 
6.3.After settling down the position and velocity of the pendulum and the velocity of 
the wheel go to zero. It is observed that by using the SFB controller the reaction wheel 
pendulum is balanced at its inverted top position with a region of attraction of ±10 
degrees. The pendulum is made stable in the vertical upward position and the wheel 
velocity goes to zero as shown in figure 6.2 The control action is shown in figure 6.3 

6.4 Fuzzy Logic Controller for Reaction Wheel Pen- 
dulum 

A fuzzy logic controller has been designed to stabilize the pendulum at its top unstable 
equilibrium point. The disk position q 2 (t) can be arbitrary, since disk position is of no 
consequence during the control. For reaching equilibrium point it is required to control 
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the position of pendulum qi(t), velocity of pendulum <ji(t) and velocity of wheel q 2 (t) . 
The disk position is left unspecified. 

The fuzzy variables are defined as follows 


e 9 i(f) = 0 

(6.15a) 

e q i(t) = 0 - g'i(t) 

(6.15b) 

e q2 (t) = 0 - q 2 {t) 

(6.15c) 


where, e q i(t) is the error in position of pendulum, e gl (t) is the error in velocity of 
pendulum and e q2 (t) is the error in velocity of wheel. 

6.4.1 Fuzzification 

For fuzzification, triangular membership functions with 9 region are chosen for the fuzzy 
variables defined in equation (6.15). Total 729 ( 9 3 ) rules are generated [12]. Maximum 
numbers of rules are fired for a given value of e q \ (t), e q \(t) and e q2 (t) are 8. 


6.4.2 Inference Mechanism 

As described in section 3.6.3, individual rule based inference mechanism is used. In this 
problem, following inference mechanism is chosen. 

Minimum of (//(gi), /i(g‘i), /r(c? 2 )) is taken as the truth value of n th rule to be fired. 
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6.4.3 Defuzzification Module (DM) 


It is described in section 2.7.3. In this problem, weighted average method [13] has been 
used for defuzzification. According to this, the control action is computed as, 


L*= 1 tk 


(6.16) 


Where the terms appeared in equation (6.16) have been described in section 2.7.3. 
Finally this value has to be denormalized so that the point wise value of the control 
output is converted to its physical value. 

Genetic Algorithm for multi objective optimization (Nsga2 code) is used as a 
optimizing tool for the fuzzy logic controller [22] . The following parameters of FLC are 
optimized using GA. 


1. Membership function for e q \ (t). 

2. Membership function for e q \ (t). 

3. Membership function for e q2 (t). 

4. Number of rules. 

Following cost functions are opted to minimize the torque and error during the simula- 
tion. 


t 


minimize ^ e 2 2 (f) 

1=0 

+ 

(6.17) 

L 

minimize r 2 (t) 

2=0 

(6.18) 
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Figure 6.4: Trajectory Tracking of Velocity and Position of Pendulum by Using FL 
Controller 



Figure 6.5: Control Torque for FL Controller with GA Optimization 
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6.4.4 Simulation Results 


The fuzzy logic controller (FLC) has been implemented for reaction wheel pendulum 
which is a nonlinear dynamical system. The simulation results are shown in the figures 
6.4 and 6.5. After settling down the position and velocity of the pendulum and the 
velocity of the wheel go to zero. It is observed that by using the FL controller the reac- 
tion wheel pendulum is balanced at its inverted top position with a region of attraction 
of ±10 degrees. The settling time is less than that of state feed back controller. After 
settling down the position and velocity of the pendulum and the velocity of the wheel 
go to zero. 


controller 

Settling time 

Control force 


pendulun 

position 

i pendulun 
velocity 

l wheel 
velocity 


SFB controller 

3sec. 

3sec 

4sec 

0.01266 

Fuzzy controller 

1.3sec. 

1.3sec 

1.3sec 

0.02815 


Table 6.1: Control of Reaction wheel Pendulum: Comparison 


Table 6.1 shows a comparative performance between state feedback controller and 
fuzzy logic controller (FLC) in terms of settling time and control input. It seen from 
this table that for state feedback (SFB) controller, settling times of pendulum position, 
pendulum velocity and wheel velocity are 3 seconds, 3 seconds and 4 seconds respec- 
tively. For fuzzy logic controller (FLC) , settling time of position of pendulum, velocity 
of pendulum and velocity of wheel are 1.3 seconds, 1.3 seconds and 1.3 seconds respec- 
tively. The root mean square (RMS) value of control torque is 0.01266 Nm for SFB 
controller and 0.012815 Nm for FL controller. 
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6.5 Summary 


The system dynamics of reaction wheel pendulum is derived using Euler-Lagrange for- 
mulation. It is a nonlinear dynamical system. A classical controller based on lin- 
earization technique and a fuzzy logic controller have been implemented to balance the 
reaction wheel pendulum in the vertical upward position. Only balancing control algo- 
rithm is implemented. The balancing is obtained with a convergence region of ±10 0 
(degrees) from vertical upward position. Fuzzy logic controller optimized with genetic 
algorithm, gives a very good performance compared with SFB controller in terms of 
settling time for same resign of attraction. 
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Chapter 7 


Conclusion and Future Work 

7.1 Conclusion 

The real world system is highly complex and nonlinear in dynamics. Studies in non- 
linear systems thus fascinated many researchers. It is difficult to obtain the exact 
dynamics of the nonlinear system as well as to design a suitable controller once the 
exact dynamical model is known. Unmodeled dynamics, parameter uncertainities are 
some of the characteristics which make the controller design more complex for nonlinear 
systems. 

Concept of fuzzy set has been proposed by Zadeh in 1965 [23] to incorporate the 
human nature of understanding. Using the concept of fuzzy set, a fuzzy logic controller 
(FLC) can be designed based on some available informations about the system to be 
controlled. No exact mathematical model is required to design the control action of 
an FLC whereas the conventional control schemes depend on the exact dynamics of 
a system. Furthermore, FLC has a capability to tolerate uncertainties in the system 
which makes it more robust compared to conventional control schemes. Thus it is 
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worthwhile to implement FLC for nonlinear systems along-with conventional control 
strategies to see the advantages of it. 

In this work, a number of nonlinear systems have been studied and various control 
techniques including FLC have been implemented on those systems. The nonlinear 
systems studied in this work include single link manipulator, reaction wheel pendulum 
and steam generator. The performance of FLC has been compared with those of con- 
ventional controllers for the above mentioned nonlinear systems. The fuzzy parameters 
have been optimized with genetic algorithm. In trajectory tracking of single link ma- 
nipulator the performance of FLC is better than PD controller and as compatible with 
controller using feedback linearization. In balancing control of a reaction wheel pendu- 
lum, the performance of fuzzy logic controller is better than state feedback controller 
in terms of settling time. 

The U tube steam generator has a nonlinear dynamics. It is required to control the 
water level of U tube steam generator for maintaining a dis-interrupted power supply. 

A fuzzy logic controller has been designed for water level control of U tube steam 
generator in a nuclear power plant using the linearized model of the plant. Its per- 
formance is compared with state feedback controller designed for the same. It is seen 
that the performance of state feedback controller is better than fuzzy logic controller 
in terms of settling time. 

For U tube steam generator, a dynamic model is necessary to analyze the perfor- 
mance of the controller. The design model should be simple enough to work with but 
must retain the essential features of the process. Detailed modeling of a U tube steam 
generator has been carried out based on the technical data provided by BARC. Some 
approximations have been made during modeling of the process. Using this model an 
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FLC, a PI and a PD controller is designed. It is seen that FLC gives better performance 
in terms of steady state error and settling time. 

The above discussion concludes that FLC can be chosen as a better alternative for 
designing the controller of a nonlinear dynamical system when the exact dynamics is not 
known. Although it gives very good result there is no specific procedure for choosing 
the rule, the membership function etc. If this drawback can be taken into account, FLC 
can be widely used for many industrial applications in near future. 

7.2 Scope for Future Work 

• The proposed fuzzy logic controller is to be tested on an actual steam generator, 
which would be my endeavor once I resume my responsibility at BARC. 

• A new FLC based on Lyapunov stability concept has been reported in the litera- 
ture. It can be implemented on the boiler model described here and a performance 
comparison can be made. 

• Many approximations have been made while modeling various components of the 
boiler. It would be interesting if these uncertain and unknown dynamics can be 
modeled using neural networks. 

• While using a linear model for boiler in chapter 3, robustness aspect has not been 
taken into consideration. Since, linear approximation may not be quite close to 
the original system, it would be worthwhile to analyze the robust stability of 
proposed controllers. 
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Appendix 


Linear Quadratic Regulator 
Consider a linear system modeled by 

x = Ax + Bu , a:(0) = xq 

with a performance index 


J = 



+ u t Ru ) dt 


where Q = Q T > 0 and R — R T > 0 

Find u = —kx such that J is minimized. 

If u is an optimum controller, then x = (A — Bk)x is also asymptotically stable. 
=$■ If V = x T Px then along the trajectory of the closed loop system is negative 
definite. Can one design an optimal control law for which ^ < —(x T Qx + u T Ru)7 This 
is proved in the following theorem. 

Theorem 

If u* = —kx such that 


min u 



+ x T Qx + u t Ru 


= 0 


for some V = x T Px, then the controller is optimal. 
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Proof : 


dv 

dt 


+ x r Qx + u* t Ru* = 0 


dv 

dt 


= -x T Qx - u* t Ru* 


Integrating 


POO 

V (rr(oo)) — V (x(0)) = — / ( x T Qx + u* T Ru *) 

Jo 


dt 


Since the closed loop system is stable, 


x(oo) 
V(x(0)) 
=r- J(u*) 


roo 

= x%Px o— / (x T Qx + u* T Ru*) dt 

Jo 


XqPxq 


To show that such a controller is optimal, here is a proof by contradiction. Let’s assume 
that 

min u + x t Qx + u T Rt?J = 0 

at u = u* , but u* is not optimal. =4> There exists some u such that J(u) < J(u*) 


dv 

dt 


+ x r Qx + u t Ru > 0 


dv 

dt 


> —x J Qx — u t Ru 
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Integrating 


V (x(0)) < / {x T Qx + u T Ru) dt 

Jo 

t T 

J(u*) < J(u) 


which is a contradiction. 

What is the optimal control law? 


f(u) = — + x T Qx + u t Ru 
at 


At the minimum value where u = u* 


g/M 

du 


— 0 and 


S 2 f(u) 

du 2 


>0 


Thus £ (| + x r Qx + vFRu) |„ u . = <F 


(2 x T Px + x t Qx + u t Ru ) 

d 

= (2 x t PAx + 2 x t PBu + x t Qx + u T Ru ) 

= 2 x t PB + 2 u t R = 0 T 
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This implies 


uFR = —x t PB 
u T = — x t PBR ~ 1 
u = -R- l B T Px=^k = R~ x B T P 

Please note that 


h{jt +xTQx+uTRv ) 

= (2 x t PB + 2u t R ) 

ou v ' 

= 2i?> 0 


The optimal closed loop system has the form 

a; = (A - BR~ l B T P)x 

However the matrix P is unknown which must satisfy the identity given in Theorem 
7.2. Thus follows the derivation for P. 


— +x t Qx + u* t Ru* = 0 

dt u=u* 

2 x t PAx + 2 x t PBu* + x T Qx + u* t Ru* = 0 
x T (A T P + PA)x - 2 x t PBR- 1 B t Px + ar T Qrr + x t PBR~ x B t Px = 0 

x T ( J 4 T P+PA + g-PSi?- 1 5 T P)a: = 0 



This is referred to as the algebraic Riccati equation (ARE). 

Please note that the control law derived here for which P is computed using ARE is 
valid only when control input is unconstrained. 

linear State Feedback Control 

It is a control technique for designing a controller for a plant with linearized model. 
Consider the state space model of a system 

x = Ax + Bu (7.1) 

The control problem can be defined as, design a state feedback gain matrix K such 
that the control law given by 

u = —Kx (7.2) 

where K G R mxn is a constant matrix. Then the closed loop system is 

x = (A — BK)x (7.3) 


The poles of the closed loop system are the roots of the characteristic equation 

det{sl n -A + BK) — 0 (7.4) 

We assume that the the designer made the selection of the desired poles of the close 
loop system, and they are 
Si, s 2 , s z s n 
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The desired closed loop pole can be real and complex. If they are complex they must 
come in complex conjugate pairs. We can form the desired close loop characteristic 
polynomial, 


a c (s) = (s - si)(s - s 2 )(s - S3) (s - s„) ( 7 . 5 ) 

— s n + o: n — 1 s n_1 + + ctis + a 0 


Our goal is to select the a feed back matrix K such that 


det(sl n — A + BK ) = s n 4 - o; n _xS n 1 + + ais + c*o (7-6) 


The above problem is called the pole placement problem. We first discuss the pole 
placement problem for the single input plants. In this case, K = k £ R lxn . The 
solution to the problem is easily obtained if the pair ( A , b ) is already in the controller 
companion form. In such case we have 


A-bk = 


0 

0 


0 

0 


0 0 0 
-a 0 — h —ai — &2 -a 2 — &3 


2 k n —i On— 1 k r 


Hence the desired gains are 
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kl = do — ao 


k2 = a\ — a\ 

kn — CX . n — i CL n — i 

If the pair (A, b ) is not controller companion form, we first transform into companion 
form and then compute the gain vector k such that 

det(sl n — A + bk) = s n + a n _xs n_1 + + axs + ao (7.7) 

thus, 

k—[cx o o 0 , Q:i Ui) cx n — i n n _x] 

Then, 

k = kT, 

Where T is the transformation that bring the pair (A, b) into controllable companion 
form. It can be represent the above formula 


kT = [of 0 - ao, ax - Ox, a n _x - a n _x][gx, q t A qiA n X ] T (7.8) 

= qi(cxoI n + cxiA + .... + cx n -iA n 1 ) — q\(a,oI n + a^A + .... + a n _i A n 1 ) 

By the Cayley-Hamilton theorem, we have 
A = — (cLoI n + Ox A + .... + a n _xA n_1 ) 
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Hence. 


k = qia c (A) 

The above expression for the gain row vector was proposed by Ackerman in 1972. 

Feedback Linearization 

The central idea of feedback linearization is to transform a nonlinear system dynamics 
into a linear one using (dynamic) feedback. The advantage is evident that once the 
system is transformed into linear one, all linear control techniques can be applied. The 
theory of linearization is developed from Differential Geometry [24]. There are basically 
two kinds of approaches for linearization, namely 

1. Input-Output linearization. 

2. Input-State linearization. 

Input- Output Linearization 

Consider a single input single output (SISO) nonlinear system 


x = f(x) + g{x)u 

(7.9) 

il 

(7.10) 


See in the above equation (7.9) output is related to input through state. The objective 
is to bring out direct linear relationship between input and output 

y m = v (7.11) 
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where v = 7 (x,u), 

m E R indicates differentiation by m times and m < n and m is known the relative 
degree of the system(7.9). 

If in the case m < n, then the transfered system (7.11) will have unobservable 
dynamics of order n — m. MIMO systems are also dealt in the same way with bit 
more complexity. This technique has been successfully implemented in several practical 
applications, such as flight control [25] and the control of rigid robots by so-called 
computed torque method [26]. 

Input-State Linearizion 

In this approach, the objective of linearizion is achieved in two steps, State transforma- 
tion and Input transformation. Moreover in this approach we are not concerned with 
output. The formula definition of the problem (in SISO case) is given in [27]. Consider 
the system 

x = f{x) + g(x)u (7.12) 

the objective is to find out a state and input transformations 


v = q(x) + s(x)u 

(7.13) 

"77 

II 

(7.14) 


so that the original system is transformed into linear one. 

z = Az + Bv (7.15) 
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In fact in the equation (7.11), relative degree, m = n, order of the system (7.9) then 
Input-Output linearization leads to Input-State Linearization with state transformed 
from x — *■ y. 

Control design using feedback linearization involves partial differential equations 
and solving these equations is really difficult particularly in the case of MIMO systems. 
Still the topic is under study and the research is underway to overcome the problem of 
solving partial differential equations. 

Under-actuated Mechanical Systems 

Under actuated systems are those systems of the form 

Q = /(?.«) +G(q)u (7.16) 


where, 

q is the state vector of independent generalized coordinates, (also called as vector of 
configuration variables). 

f(.) is the vector field representing the dynamics of the systems. 
q is the generalized velocity vector. 

G is the input matrix. 

and u is a vector of generalized force inputs. 

The dimension of q is defined as the degrees of freedom of the system . In general (rather 
literally) a system is said to be under actuated if the external generalized forces are 
not able to command instantaneous accelerations in all directions in the configuration 
space, i.e. rank(G) < dim(q). 
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(7.18) 

(7.19) 

(7.20) 

(7.21) 

(7.22) 

(7.23) 

(7.24) 


From equation (4.10) and (7.24), it proves the formula used for calculation the 
homogeneous density. The parameters used from boiling section to U tube resign are 
V— Volume of steam generator. 
m,f = Mass of flow(Kg). 
m s = Mass of steam flow (Kg). 
m w = Mass of water flow (Kg). 

V 3 = Volume occupied by steam (Kg). 
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V w = Volume occupied by water meter 3 . 


Terms related to modeling of the process 

Saturation temperature: The temperature at which a liquid substance boils under 
a given pressure. 

Latent heat of vaporization: The heat that changes the physical state of a sub- 

/ 

stance from a liquid to a vapor, or from a vapor to a liquid; no temperature change is 
shown by a thermometer during the conversion process. It is a function of pressure and 
temperature. 

Specific heat: The quantity of heat require to raise the temperature of a substance 
of unit mass to unit degree under a given pressure and temperature. 

Saturated steam: Steam at the temperature and pressure at which evaporation oc- 
curs. 

unsaturated steam: when both the water and liquid are present steam and water 
mixture, then the steam is called unsaturated steam. Saturated water: Water at its 
boiling point. 

Steam quality: The percent by weight of vapor in a steam and water mixture. 
Density : It is the ratio of mass to volume of a substance. 

Homogeneous density : When the density is calculated for the combination of two 
substance it is called homogeneous density, as described in section 4.5.1, the homoge- 
neous density for water and steam is calculated. 
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